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A major problem confronting the understanding of tidal evolution of close solar-type binaries is the inefficiency of the 
turbulent convection. The value of the effective viscosity estimated, in the framework of the mixing length theory (MLT), 
implies circularization timescales which are almost two orders of magnitude longer than observed. Moreover, the reduction 
of the effective viscosity due to the fast time- variation of the tidal shear in short period binaries, increases the discrepancy 
to about three orders of magnitude. This state of affairs has motivated suggestions that tidal orbital evolution, notably 
circularization occurs mainly during the pre-main-sequence phase. However, observational data accumulated over the 
recent decades imply that circularization does occur during the the main-sequence phase (Mazeh 2008). 
In this work, we examine the possibility that the apparent inefficiency of turbulent convection is merely a shortcoming of 
MLT approach. Indeed, a recent 3D numerical simulation (Penev et al. 2007), suggests that the true convective viscosity is 
probably larger than the MLT value and that the reduction due to the time- variation of the shear is not drastic. We employ 
a model for stellar turbulent convection (Canuto, Goldman & Mazzitelli 1996) to evaluate the effective viscosity both for 
a steady for and time dependent tidal shear. 

The model is physically based, self-consistent, and accounts for the full spectrum of the turbulent eddies. It has been found 
advantageous, compared to the MLT, in many applications. We use an analytic approximation to the turbulent spectrum to 
obtain the reduction of the efficiency due to the time- variation of the tide. The results are: 

(i) an enhanced effective viscosity ( by a factor of ~ 4.5) and more importantly 

(ii) only a mild reduction due to the time-variation of the tidal shear. 

Overall, for binaries with orbital period of 15 days the discrepancy is "only" a factor of ~ 30 down from a factor of 
~ 1000. These encouraging results should motivate an investigation of rigorous non-analytic solutions. 
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1 Introduction 

About two decades ago there has been a revival of interest in 
the tidal orbital evolution of close binary systems. This was 
mainly due to the availability of large samples of spectro- 
scopic binaries, made possible by the (then) new generation 
of stellar speedometers in full operation at several obser- 
vatories (Latham et al. 1988; Mathieu and Mazeh 1988); 
Latham et al. 1992a; Latham et al. 1992b; Mathieu 1992; 
Mathieu et al. 1992). More recently, the interest in the sub- 
ject has been renewed by Latham et al. (2002), Mathieu et 
al. (2003), Mathieu (2005), Mazeh (2008), Meibom (2005), 
Meibom and Mathieu (2004, 2005), Meibom et al. (2006), 
Meibom et al. (2007). 

From the theoretical point of view, the central issue is 
the nature of the mechanism responsible for the dissipation 
of the energy in the tide and for the transfer of angular mo- 
mentum between the binary orbit and the stellar rotations. 
A successful model should predict the dependence of ob- 
served parameters the binary orbital parameters. 

For these binaries, it is commonly believed that the dom- 
inant physical process controlling energy dissipation is the 
turbulent convection in the stellar envelop. In this equilib- 
rium tide model (Alexander 1973; Lecar et. al., 1976; Zahn 
1977; Hut 1981, 1982; Eggleton et al. 1998), the tidal force 




induces a sheared velocity field in each of the stars. The 
turbulent convection interacts with the shear drawing from 
it energy. Technically, the effect of the turbulence is repre- 
sented by an enhanced effective diffusivity- the turbulent 
viscosity which is orders of magnitude larger than the mi- 
croscopic viscosity. Within this framework, the circulariza- 
tion timescale for late-type stars with convective envelopes 
is proportional to the orbital period of the binary to the 
power 16/3 (Zahn 1977). 

Zahn (1989) called attention to the fact that the tidal 
effective viscosity is reduced because the orbital period is 
shorter than the typical timescale of the turbulent convec- 
tion largest eddies. Thus, the effect of the large eddies is av- 
eraged out and only smaller eddies, with shorter timescales 
can contribute to the viscosity. This implies a reduction of 
the viscosity and thus implies an increase of the circular- 
ization timescales. According to Zahn (1989) the reduction 
is linear in the period thus implying that the circularization 
timescale is proportional to the orbital period to the power 
13/3. As a result, Zahn and Bouchet (1989) argued that most 
of the orbital circularization should occur during the pre- 
main-sequence phase. 

Goldman & Mazeh (1991) concluded, in analogy with 
kinetic microscopic viscosity, that the reduction should be 
quadratic in the period implying that the circularization 
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timescale is proportional to the orbital period to the power 
10/3. Such a quadratic reduction has been proposed by Gol- 
dreich & Keeley (1977) in the context of the damping of 
solar pulsations by turbulence convection. The observations 
available at the time favored an index of 10/3. This led Gold- 
man& Mazeh (1992) to suggest that the model is essentially 
correct and the excessively long timescales reflect the MLT 
shortcomings, notably not representing the true topology of 
the turbulent convection and not accounting for the entire 
turbulent spatial energy spectrum. Moreover, since the MLT 
is basically a descriptive framework, the expressions used 
for the turbulent viscosity and for the turbulent timescales 
are based on dimensional analysis and the numerical values 
of the various parameters are thus uncertain. 

The reduction of the turbulent convective viscosity has 
been readdressed by Goodman and Oh (1997). They ob- 
tained that the reduction increases gradually as the orbital 
period gets shorter; it is quadratic in the limit of very short 
orbital periods. Using the MLT, they found that the ensuing 
circularization timescales for an orbital period in the range 
of 15 - 10 days are longer by an additional factor of 10. 
Overall, for a 15 day binary period, the predicted timescale 
is too long by a factor close to 1000. 

Recent observations (Meibom & Mathieu 2004; 2005; 
Meibom et al. 2006; Meibom et al. 2007; Mazeh 2008) re- 
veal a more complex picture than that of the past. Tidal 
evolution does take place during the main sequence phase. 
However, the functional dependence of the tidal timescale 
is not clear; the observational uncertainties do not allow a 
determination of the index of the power-law, if indeed there 
is a description by a single power-law. 

An interesting paper by Penev et al. ( 2007) reports the 
results of a 3D numerical simulation suggest that the true 
convective viscosity is probably larger than that indicated 
by the MLT for the relevant orbital periods. 

2 Motivation for Revisiting the Problem 

The new observational data mentioned in the last paragraph, 
seem to imply that circularization does occur during the 
main sequence phase. Since, at least qualitatively, turbulent 
convection is the natural candidate for supplying the effec- 
tive viscosity, we wish to reconsider the problem in a theo- 
retical framework more consistent than the MLT. Support in 
this direction is given by Penev et al. ( 2007). 

We use a model for turbulent convection of Canuto, 
Goldman, & Mazzitelli (1996) (CGM96) which is based on 
previous turbulence models (Canuto et al 1984; Canuto and 
Goldman 1985; Canuto et al. 1987; Canuto et al. 1988). The 
CGM96 model has been applied to a large body (over 130 
citations in ADS) of astrophysical problems, spanning a va- 
riety of aspects of stellar convection: stellar structure, stel- 
lar evolution, age of globular clusters, helio seismology, and 
blue-edge of DBV pulsating white dwarfs . 

The CGM96 model is a self-consistent model for turbu- 
lent convection which accounts for the full spectrum of the 



turbulent convective eddies. It constitutes a significant im- 
provement over the mixing length theory (MLT) which is a 
one eddy approximation to the full energy spectrum of the 
turbulence. The model is local and homogeneous and does 
not account explicitly for compressibility and anisotropy; 
these could however be addressed approximately. Techni- 
cally, it consists of a set of coupled ordinary differential 
equations which are easily solved numerically. 

The turbulence energy spectrum is determined by the 
rate of energy input from the source (buoyancy) into the 
turbulence. This rate in turn depends both on the source pa- 
rameters and on the turbulence itself. The physical basis of 
the model is probably the reason for its success in predicting 
observational results, as mentioned above. 

The model complements numerical simulations by pro- 
viding physical understanding and an easy tool to test for 
the imprints of the various parameters of the problem. It can 
be also used a sub-grid model for Large Eddy Simulations 
(LES). In the context of the present work we note that 

1 . Since CGM96 includes the full turbulence spectrum, the 
turbulent viscosity is expected be larger than the MLT 
value. 

2. Since CGM96 has definite prescriptions for the turbu- 
lent viscosity and for the convection timescales - the 
ambiguity of the MLT prescription could be avoided. 

In what follows we report on the results of this investiga- 
tion. A more detailed account will be presented in Goldman 
(2008). 

3 The Effective Turbulent Viscosity in the 
Presence of Time- Varying Shear 

Let us consider stationary turbulence that interacts with an 
external shear. We focus on the case when the rate of en- 
ergy input from the shear to the turbulence is small com- 
pared to the energy rate fed into the turbulence by its source 
(buoyancy in the present case) and does not modify the tur- 
bulence. 

The turbulence is defined by F(k), the turbulent veloc- 
ity spectrum 

/>OC 

< v(t) -v(t)) >= / F(k)dk (1) 

J ho 

which is related to turbulent energy (per unit mass) spec- 
trum, E{k), by F(k) = 2E(k). 

The wavenumber k corresponds to the largest scale in 
the turbulence energy spectrum with F(k ) = 0. The spec- 
tral function increases with k and then decreases and ap- 
proaches a Kolmogorov spectrum, F(k) oc A; -5 / 3 . For large 
enough k, microscopic dissipation truncates the spectrum. 

The quantity that controls the interaction between the 
turbulence and the shear is the turbulent diffusivity Dt- It 
relates the turbulent stress tensor (caused by the external 
shear) with the shear tensor. 

<v i {t)v j {t)>=D T S ij {t) (2) 
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The angular bracket denotes ensemble average; the turbu- 
lence is stationary, hence the result is independent of t. The 
turbulent diffusivity Dt, controls also turbulent transport 
processes. 

The turbulent diffusivity Dt has the same dimensions 
as the turbulent eddy viscosity vt- However, the latter con- 
trols the non-linear interactions of the eddies, while the for- 
mer controls the interactions with an external shear. Given a 
specific turbulence model, these two quantities are defined 
and the relation between them could be established. In a 
phenomenological approach, like the MLT, the two are con- 
sidered to be the same and taken equal (up to some dimen- 
sionless coefficient) to the ratio of the eddy size to the tur- 
bulent velocity on this scale. In the CGM96 model these are 
distinct concepts. 

For the problem at hand Eq. (O takes the form 

<V r (t)v <f> (t)>=D T Sr^(t) (3) 

Using the shearing coordinates approach (Rogallo 198 1) 
we obtain 

< v r {t) Vrl> {t) >= / < v r (t)v r (t') > S r4> (t')dt' (4) 




Fig. 1 Solidline: the efficiency factor f](x). Thinlines: 
x , x~ 2 , representing linear and quadratic reductions, 
thatt 



3.1 Analytic Approximation of the Efficiency Factor 

The efficiency factor can be evaluated once the turbulence 
spectrum and the eddy decorrelation time scale are known. 

If the turbulence spectrum of CGM96 is approximated 
by a Kolmogorov spectrum (even for the largest eddies), it 
is possible to obtain an analytic approximation: 



The 2-times ensemble average is related to the turbu- 
lence energy spectrum by (Monin & Iaglom 1975; Voelk et 
al. 1980) 



r)(x) 



20 



-ln(l + x 2 ) - - + 



x 



^ ( i-^4M))(ii) 



X 



< v{t) ■ v{t')) >= / F(k)e dk 

Jko 

For the quantity of interest in Eq. © we get 



(5) 



< v r {t)v r (t') >= / a(k)F(k)e dk (6) 
Jk 

In the expression above, r c (fc) denotes the eddy decorrela- 
tion timescale and < a(k) < 1 , the anisotropy factor of 
the eddies defined by 



< v r (t)v r (t) >= / a(k)F(k)dk. 
Jk 



(7) 



For strictly isotropic turbulence a(k) = |. 
For a pure sinusoidal shear S(t) = Se luJt we get an 
expression for the diffusivity: 

D T (w) = I" a{k)F{k)e~^ {k) ^ jdfc (8) 

Jk 1 + [uJT c {k)} 

For slow varying shear (uiT c (k) << 1, for all k), the tur- 
bulent diffusivity contributed by the full turbulent spectrum 



The efficiency function is shown in figure Q] together 
with efficiencies that decrease as , and x~ 2 , respec- 
tively, normalized to unity at x=0.5. As seen, the decrease 
with x is quite moderate. The qualitative shape is similar to 
that in Goodman and Oh (1997). 

4 Application to the Solar Model 

We have applied the approach outlined above to a binary 
consisting of two solar mass stars, and used the turbulence 
model of CGM96 to compute the relevant quantities. The 
details will be published elsewhere (Goldman 2008). The 
main results are presented below. 

4.1 Tidal Timescales for a Stationary Shear 

The basic tidal timescale T c the circularization time, is in- 
versely proportional to a weighted average, over the convec- 
tive zone, of the product pDx'- 



(12) 



Dt(0) 



a(k)F(k)T c (k)dk 



(9) 



When the shear is varying faster than the decorrelation 
timescale of the largest eddies, r , the diffusivity is reduced 
and can be expressed as 



< pD T >= / pD T x a dx 

Jo.7 

with x = r/ Rq . 

For a stationary shear 

< pD T {0) >= [ pD T (0)x 8 dx 

Jo.7 



(13) 



Dt(lo) = T)(ujto)D t (0) 
where i](ujT n ) < 1 is the efficiency factor. 



(10) 



The numerical solution of the solar model of CGM96, 
yields 

< pD T (0) >= 2.4 x 10 11 grcm 2 s- 1 (14) 
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This is is a factor of ~ 4.5 larger than the corresponding 
value computed within the MLT framework. For a binary 
consisting of solar mass members, with orbital period of 15 
days the circularization timescale is T c (0) = 1.7 x 10 11 yr. 
So, the discrepancy with the observational circularization 
timescale is alleviated by the same factor and is now "only" 
a factor of ~ 12 too large. 

The contribution to < pDx (0) > peaks in the middle of 
the convective zone, reflecting the facts that the density is a 
decreasing function of the radius and the diffusivity Dt{0) 
is an increasing function. 

The increase of < pDt(0) > in comparison to the MLT, 
is not as large as might have been anticipated. According to 
CGM96, taking into account the full turbulent spectrum en- 
hances the convective flux for a given super adiabatic gradi- 
ent, by an order of magnitude, compared to the MLT. Since 
the value of the convective flux is imposed by the radiative 
core, the super adiabatic gradient is smaller, than in the MLT 
models. This compensates partially, the increase stemming 
from the inclusions of the full turbulence spectrum. 

4.2 Tidal Timescales for a Dynamic Shear 

We have computed at each radius the local eddy decorre- 
lation time scale To, and from it the local efficiency factor 
which was then used to evaluate < pDr > from Eq. [T2] 

We obtained that for a 15 day binary, the weighted ef- 
ficiency is ~ 0.4. This very mild reduction leads to a tidal 
time scale which is too long by a factor of ~ 30, compared 
to a factor of ~ 1000, according to the MLT. 

5 Discussion 

The results indicate that the use of a physically based model 
for turbulence alleviates considerably the discrepancy be- 
tween the observational and predicted tidal time scales. The 
approximation, used here, of a Kolmogorov spectrum even 
for the largest eddies, underestimates the diffusivity. There- 
fore, we intend to perform a more rigorous, non-analytical 
computation. 
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